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Abstract 

Based on the Kaplan-Nelson Lagrangian, we investigate kaon condensation in dense 
neutron star matter allowing for the explicit presence of hyperons. Using various models 
we find that the condensate threshold is sensitive to the behavior of the scalar density; the 
more rapidly it increases with baryon density, the lower is the threshold for condensation. 
The presence of hyperons, particularly the E~ , shifts the threshold for K~ condensation 
to a higher density. In the mean field approach, with hyperons, the condensate amplitude 
grows sufficiently rapidly that the nucleon effective mass vanishes at a finite density and 
a satisfactory treatment of the thermodynamics cannot be achieved. Thus, calculations of 
kaon-baryon interactions beyond the mean field level appear to be necessary. 



The idea that, above some critical density, the ground state of baryonic matter might 
contain a Bose-Einstein condensate of negatively charged kaons is due to Kaplan and Nelson 
12). Subsequently, the formulation, in terms of chiral perturbation theory, was discussed 

and the astrophysical consequences have been explored in ref. f|]. Physically, the 
strong attraction between K~ mesons and nucleons increases with density and lowers the 
energy of the zero-momentum state. A condensate forms when this energy becomes equal 
to the kaon chemical potential, /x. In catalyzed (neutrino-free) dense neutron star matter, 
/i is related to the electron and nucleon chemical potentials by fi = \x e = p-n — Hp due to 
chemical equilibrium in the reactions n <-> p + e~ + v e and n «-> p + K~ . The density at 
which this takes place is model and parameter dependent, but is typically ~ 4no, where no 
denotes equilibrium nuclear matter density. Since this may be less than the central density 
in neutron stars, a K~ condensate is expected to be present in the core region. Apart from 
the softening effect on the equation of state, which lowers the maximum mass, the proton 
abundance is dramatically increased resulting in a nucleon, rather than a neutron, star. 

However, many calculations, e.g. || q], of dense matter indicate that hyperons, starting 
with the S~ and A, begin to appear at densities ~ (2 — 3)no> a little lower than the kaon 
threshold mentioned above. Once a significant number of negatively charged hyperons are 
present, less electrons are required for charge neutrality, so the electron chemical potential, 
He, begins to decrease with density. Since [i e = \x this would be expected to inhibit kaon 
condensation, delaying it to a higher density. However, since overall charge neutrality must 
be maintained in the presence of many strongly interacting particles, the issue of whether 
or not kaon condensation would occur in matter containing hyperons is uncertain. The only 
previous work JjJ involving both hyperons and kaons was concerned with p-wave interactions 
which led to quasinucleons through the mixing with hyperons. 

Our objective here is to explore the influence of the explicit presence of hyperons on 
kaon condensation in dense neutron star matter. In our model calculations, we employ the 
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Kaplan-Nelson Lagrangian for the kaon-baryon interactions and a relativistic field theoret- 
ical approach for the baryon-baryon interactions. Use of the relativistic approach for the 
baryons allows for a distinction between the scalar and vector densities, which is necessary 
for an adequate treatment of the kaon-baryon interactions. Ignoring this distinction results 
in a lower threshold density for kaon condensation. 

We begin with the Kaplan-Nelson SU(3) x SU(3) chiral Lagrangian for the kaons and 
the s-wave kaon-baryon interactions. Specifically 

C K = \f 2 Trd^Ud> M U+ CTim q (U + U ] -2) + %TiB'f\y ix , B] + aiTrB(fm,f +h.c.)B 
+a 2 TrBB(£m q £ + /i.e.) + a 3 {Tr BB} {Tr(m q U + h.c.)} . (1) 

Here U is the non-linear field involving the pseudoscalar meson octet from which we retain 
only the contributions, £ 2 = U and B represents the baryon octet - nucleons plus 
hyperons. The quark mass matrix m q = diag(0, 0, m s ). For the mesonic vector current, Vn, 
only the time component survives in an infinite system with V = K&do£ + {dog). The 
pion decay constant / = 93 MeV and C, a\,a^ and 03 are constants. After some algebra, 
the relevant part of Ck takes the form 

V X J { COS 2 ±X b 



-[ m K + 7TH^l[( ai + ° 2 )( 1 + y bQb) + (01 - a 2 )(qB - Yb) + 4a 3 

^ J B 



BB 



+||(ai + a 2 )(2AA + V3[t A + AS ])) , (2) 

OJ J COS 9 X 



m. s , s,„ T . . /^n . T^nnx\ K + K 

cos 2 \x. 

where x 2 = 2K + K~ / f 2 , qs and Yb are the baryon charge and hypercharge, respectively, 
and the kaon mass is given by m 2 K = 2Cm s / f 2 . We have not included in Eq. (||) terms which 
simply give a constant shift to the baryon masses; they indicate that aim s = — 67 MeV 
and a 2 vn s = 134 MeV using the hyperon-nucleon mass differences. The remaining constant 
a^m s is not accurately known and we shall use values in the range —134 to —310 MeV 
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corresponding to to 20% strangeness content for the proton. Some guidance is provided 
by recent lattice gauge simulations || which find that the strange quark condensate in the 
nucleon is large, i.e. (JV|ss|JV) = 1.16 =t 0.54. Prom the relation m s (~ss) p = — 2{a 2 + a^)m s 
and using m s = 150 MeV, we obtain a^m s = —(220 ± 40) MeV, which is in the middle of 
our range of values. 

In the baryon sector, we employ a relativistic field theory model in which baryons 
interact via the exchange of a-, p- and w-mesons. In the case that only nucleons are 
considered, B = n,p, this is the well-known Walecka model || and we allow for possible 
"non-linear" <r 3 and <r 4 terms. The Lagrangian takes the form 



C B =J2 B (i^d^ - g uB -fw» ~ gpB-fTop -t-M B + g aB a)B - \F^ V + ±m, 



B 



-±B^ • B"" + \m 2 p h» • W + \d^u - \m 2 y - \bM{g m af - \c(g an of. (3) 

Here M B is the vacuum baryon mass, the p-meson field is denoted by b^, the quantity 
t denotes the isospin operator which acts on the baryons and the field strength tensors 
for the vector mesons are given by the usual expressions:- F^ v = d^ojy — d^u^, = 
dubv — dybu- The nucleon mass, M = 939 MeV, is included in the penultimate term so 
that b is dimensionless. 

The total hadron Lagrangian is then Ctot = £-K + We shall treat the kaons in the 
mean field approximation. For the baryons, we shall consider calculations at the mean field 
level (with the "non-linear" terms so as to obtain a reasonable compression modulus for 
nuclear matter) and at the one-loop Hartree level (without the "non-linear" terms). As 
we shall see, there is a qualitative difference between the results. We need to calculate the 
potential, $7, of the grand canonical ensemble at zero temperature. Notice first that the 
A — S° mass matrix needs to be diagonalized (for details see ref. ||) producing eigenstates 

(l + «5 2 )^ (l + 5 2 )§ yJ 
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with corresponding masses Mjj and Mg ■ . These are henceforth included in the sum over 
baryon states B, along with the n, p, rP~ . The kaon-baryon interactions can be 

included in the baryon effective masses and chemical potentials. The remaining kaon terms 
are then easily treated by writing || the time dependence of the fields = -^fOe^ 1 ^; 
thus, 9 gives the condensate amplitude. Then it is straightforward to obtain 



^ = f 2 {2m 2 K sm 2 \( 



sin 2 #) + \m 2 a a 2 + \bM(g an af + \c(g an a) 



1^,2 , ,2 l™2t2 



in 



B 



J dkk 2 (E* B -v B ) + AE{M* B ) 



(5) 



Here V is the volume, E* B = ^jk 2 + M B 2 and the effective masses, M B for Hi and H 2 are 
obtained by diagonalizing the mass matrix, as we have mentioned, while the remaining 
cases are given by 

M* B = M B - g uB a + [(ax + a 2 )(l + Y B q B ) + (ai - a 2 ){q B - Y B ) + 4a 3 ]m s sin 2 ±# . (6) 
The chemical potentials hb are given in terms of the effective chemical potentials, v B , by 
= + QujBUq + g P BhBbo - (Y B + qB)n sin 2 i6> , (7) 
where t% B is the z-component of the isospin of the baryon and the relation to the Fermi 



h 2 4- M* 2 
K FB + 1V1 B ■ 



momentum kp B is provided by u B - 

In the one loop Hartree approximation (with b = c = 0) , an additional contribution 
AE(M B ) to the energy density is introduced due to the shift in the single-particle energies 
caused by the negative energy baryon states. After removing divergences, AE(M B ) can be 



written [10] in the form 
1 

8tt 



AE(M B ) = --j 



' A ' ' M + ln m) Mb(Mb ~ Mb? ~ ln M {Mb 



m b y 



M* 

+M* B 4 ln -r + M%(M B - M* B ) - \M 2 B [M B - M* B ) 2 
M B 

+fM B (M B - M* B f - || (Ms - M* B ) 



(8) 
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Here the necessary renormalization introduces a scale parameter, \i r . For the standard 
choice |J of /j, r /M=l, the first two terms in Eq. (g) vanish. Other choices of \i r jM introduce 
explicit cr 3 and <7 4 contributions. The freedom to vary fj, r /M allows the density dependence 
of the energy to be modified, which makes it possible pi to explore equations of state 
with different stiffnesses. We refer to this approach as the modified relativistic Hartree 
approximation (MRHA). In previous work [jO] ], without hyperons, we found that while 
neutron star masses do not significantly constrain /i r /M, finite nuclei favor a value of 0.79. 

The thermodynamic quantities can be obtained from the grand potential in Eq. (||) in 
the standard way, thus the baryon number density ng = (37r 2 )~ 1 kp B , while for kaons 

n K = / 2 (/isin 2 # + 4esin 2 ±#) with e = ^T(Y B + q B )n B /(4f 2 ) • (9) 

B 

The pressure P = — Cl/V and the energy density e = —P + J2b ^b^b + M n if- The meson 
fields are obtained by extremizing S7, giving 

m l u o = 9uiB n B , m 2 p b = Y 9pbHbub and 

B B 

m 2 a a = -bMgl N a 2 - cg* N a 3 + £ g aB n s B . (10) 

B 

Here n s B denotes the baryon scalar density 

kpg 

The condensate amplitude, 9, is also found by extremizing f2. This yields the solutions 
9 = (no condensate), or, if a condensate exists, the equation 

fi 2 cos 9 + 2e/i — m 2 K — d\ — d 2 = where , 

2fdi= [(ai + a 2 )(l+Y B q B ) + (a 1 -a 2 )(qB-yB)+^a 3 \m s n s B , 

B^H U H 2 

2/ 2 d 2 sin 2 ±0 = J2 (M* B + gaAa)n B - M A n s Hi - Mxn s H2 
B=H 1 ,H 2 

+ TT5 2 ~ {Ma ~ Ms)(n ^ " ' (12) 



where we have taken g a \ = g^s ■ 

Finally we consider the leptons. The electrons and muons are governed by chemical 
potentials fj, e = fJ,u = /i. (In the cold catalyzed state, the antineutrinos would have left the 
star.) Their contributions to the total energy density and pressure are given adequately by 
the free gas expressions. The requirement of chemical equilibrium in the weak processes 
yields 

M#i = Mtf 2 = MS = Mn ; MS- = MS- = Mn + Me Mp = ME+ = Mn ~ Me • (13) 

The remaining condition is that of overall charge neutrality, namely 

1 bUb ~n K -n e -n^ = . (14) 

B 

The couplings g an and g un are determined by fitting the equilibrium empirical properties 
of nuclear matter, while the coupling g pn is obtained by fitting the nuclear symmetry energy. 
For the mean field model, we use the last two sets of values in table 2 of ref. [O]. These 
models, termed models A and B hereafter, have compression moduli of 300 and 240 MeV, 
respectively, and therefore exhibit different stiffnesses for the high density equation of state. 
The couplings employed for the MRHA model are shown in table 1. 

Large uncertainties exist in the values of the hyperon-meson coupling constants. We re- 
duce the number of parameters by making the reasonable assumption that all the hyperon- 
meson coupling constants are the same as those of the A. Defining x a = g a A/9cm, with 
analogous definitions for the u— and /^-couplings, the binding energy of the lowest A level 



in nuclear matter at saturation yields [12] 



-28 = x w g wn uJo - x a g an a , (15) 



in units of MeV. In ref. [12], on the basis of fits to hypernuclear levels and neutron star 
properties, the value x a = 0.6 was suggested, is then determined and the choice x p = x a 
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was made. We adopt this procedure, since we have examined the effect of varying x a and 
setting x p = x u and found similar qualitative behavior. 

To highlight the influence of hyperons on kaon condensation, it is useful to first consider 
the case where they are absent, i.e. nucleons-only matter. Results for this case are contained 
in table 2, and, in figures 1 and 2. The threshold density n c for condensation is determined 
from 



M H r^, p - m K 



, )[2 r L 2oinJ + (2a 2 + 4a 3 )(nJ + <)J — | = 0, (16) 

which is obtained by setting = in Eq. fll2] ) and restricting the sum over B to n and p 
only. Expressing ratio to equilibrium nuclear matter density, no, we show in table 

2 results for u c = n c /n,Q. The larger the magnitude of a^m s , the lower is the value of u c . 
In addition, u c is systematically lower for the MFT models than for the MRHA models. 
This latter feature is a consequence of the more rapid increase of the scalar densities in 
the mean field model than in the MRHA model (contrast panel (c) in fig. 1 and fig. 2). 
Due to the rapidly growing condensate amplitude (see panel (e)), the net negative charge 
carried by the kaon fields grows rapidly with density so that the leptons play only a minor 
role in maintaining charge neutrality (the partial fractions are shown in panel (a)). At high 
density, matter contains nearly equal abundances of neutrons and protons. The softening 
induced by the presence of kaons has the consequence that the maximum masses of stars are 
reduced from the case without kaon condensation |], |3| . It is worthwhile to point out that if 



the distinction between the scalar and vector densities is ignored in Eq. (16), condensation 
sets in at somewhat lower densities than those shown in table 2. 

To summarize the case of nucleons-only matter: the onset of kaon condensation is 
controlled by two factors, the first being the behavior of the scalar density in dense matter 
and the second being the magnitude of a^m s which is related to the strangeness content of 
the nucleon. Additional work is needed to further pin down these quantities. 
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We turn now to the role of hyperons in matter, results for which are shown in fig. 
3 for the mean field model B with a^m s = —222 MeV. Panel (a) of this figure shows 
the composition. One expects that the A and the E~ first appear at roughly the same 
density, because the somewhat higher mass of the E~ is compensated by the presence of 
the electron chemical potential in the equilibrium condition of the E~ (see Eq. (|l3|)). More 
massive and more positively charged particles than these appear at higher densities. Use 
of Eq. ( |T5| ) to constrain the hyperon couplings has the consequence that the first strange 
particle to appear is the negatively charged S~ hyperon. The £~ competes with the leptons 
in maintaining charge neutrality, so the lepton concentrations begin to fall. This is reflected 
in the electron chemical potential fi e (see panel (b)) which saturates at about 200 MeV and 
begins to decrease from this value with increasing density. Since fi e = fi, this has the 
consequence that kaon condensation now takes place, i.e. Eq. (|l2| ) is satisfied, at a higher 
density than in nucleons-only matter. This effect is seen for all the cases listed in table 2. 
Some insight into the role of the scalar densities may be gained by examining the threshold 
condition Eq. ( |l2|) in the presence of the S~ and A hyperons only: 



2 (2n p + n n -n E -) 2 

M + y2 ^ ~ m K 



2am s p + (2a 2 + 4a 3 )(nj; + n s n + 



m. 



+ (§(«! + a 2 ) +4a 3 )<j —| = 0. (17) 

The first two terms in this equation are smaller than in the nucleons-only case and this has 
to be compensated by the last term, which requires a higher density. The increase in the 
critical density is less in the MFT models, since the scalar densities of the various particles 
increase more rapidly than in the MRHA models. In fact, for the MRHA the central density 
of a star containing all the baryons, but no kaons, (u cen t m table 2) is less than the critical 
density u c in almost all cases, which implies that kaons will not be present. This is also the 
case for the MFT when the magnitude of a^m s is the smallest; however for larger values, 
which we favor, a kaon condensate will be present in the star. 
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Returning to fig. 3, which corresponds to the MFT with a^rris = —222 MeV, we see 
that in the presence of hyperons, the condensate amplitude increases rapidly with density 
(panel (d)); so rapidly in fact that large changes are induced in the scalar densities (panel 
(c)) and the Dirac effective masses of all the particles (panel (b)). The kaon contribution 
to the Dirac effective masses of the nucleons 

M* = M - g an a + (2a x + 2a 2 + 4a 3 )m s sin 2 ^ , 

M* = M - g an a + (2a 2 + 4a 3 )m s sin 2 \0 , (18) 

is negative and this causes the effective masses to vanish at a finite baryon density. Thus, 
within the mean field scheme adopted here and without an explicit treatment of a hadron 
to quark matter transition, a calculation of the thermodyamics to encompass the full range 
of densities in a neutron star is precluded. This indicates the need to consider several 
improvements: (i) inclusion of additional terms in the effective Lagrangian, such as the p- 
wave interactions, (ii) inclusion of loop corrections in kaon-baryon interactions, in particular 
an exact evaluation of the zero-point energy with the non-linear Kaplan-Nelson Lagrangian. 

In conclusion, we have shown that in dense matter the presence of hyperons leads to 
kaon condensation at higher densities than would otherwise be the case. The critical den- 
sities depend sensitively on the behavior of the scalar densities or, equivalently, the baryon 
Dirac effective masses in dense matter. A more rapid variation in these quantities results in 
lower thresholds for condensation. A mean field treatment of the kaon-baryon interactions 
produces rather large effects so that the nucleon effective masses vanish at a finite baryon 
density, indicating unphysical behavior. Clearly further work is needed to fully understand 
the role of the hyperon degrees of freedom on kaon condensation in dense stellar matter. 

We thank Gerry Brown for much encouragement and for valuable discussions. We are 
grateful to Norman Glendenning for much helpful correspondence. This work was supported 
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Figure Captions 

Fig. 1. Results are for nucleons-only matter in the mean field model B with a^rris = —222 
MeV. (a) Relative fractions Yi = riij \n n + n p ). (b) Nucleon Dirac effective masses, the kaon 
chemical potential \i = \i e and the scalar field a. (c) Scalar densities with n s = n s n + n p . 
(d) Pressure with and without kaon condensation, (e) Condensate amplitude 6 in degrees. 

Fig. 2. Results are for nucleons-only matter in the Hartree MRHA model with fi/M = 0.79 
azm s = —222 MeV. (a) Relative fractions Yi = m/(n n + n p ). (b) Nucleon Dirac effective 
masses, the kaon chemical potential \i = [i e and the scalar a field, (c) Scalar densities with 
n s = n s n + n p . (d) Pressure with and without kaon condensation, (e) Condensate amplitude 
9 in degrees. 

Fig. 3. Results are for matter with hyperons in the mean field model B with asm s = —222 
MeV. (a) Relative fractions Yi = ni/(J2B n B)- (b) Baryon Dirac effective masses, the 
kaon chemical potential \i = /j, e and the scalar field a. (c) Baryon scalar densities, (d) 
Condensate amplitude 6 in degrees. 
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Table 1 
MRHA coupling constants 





Without hyperons 




With hyperons 




M 


K 


r 2 


cl 


r 2 


K 


r 2 


cl 


r 2 




0.79 


354 


180.6 


317.5 


73.5 


177 


118.7 


258.1 


84.8 


0.66 


1.00 


461 


137.7 


215.0 


81.6 


455 


133.1 


210.3 


82.4 


0.66 


1.25 


264 


78.6 


178.6 


90.8 


228 


64.9 


174.0 


92.6 


0.68 



Here Cf = (giNM/mi) 2 and Kq (in MeV) is the compression modulus. The 
parameters are fitted to a nuclear matter binding energy of 16 MeV and a 
symmetry energy of 30 MeV at a density of 0.16 fm -3 . For the hyperon case 
we take x a = x p = 0.6. 

Table 2 



Critical density ratio, u c = n c /no, for kaon condensation 









Without hyp 


erons 




With hyperons 








a 3 m s 


-134 


-222 

u c 


-310 


-134 


-222 


-310 


a) 
w ccnt 


MFT 


A 




4.14 


3.14 


2.49 


8.97 


4.20 


2.74 


6.72 




B 




4.15 


3.15 


2.49 


9.46 


4.22 


2.73 


7.66 




Hr/M = 


0.79 


5.22 


4.54 


3.89 


15.3 


9.87 


7.32 


7.81 


MRHA 


Hr/M = 


1.00 


5.12 


4.33 


3.61 


13.9 


8.64 


6.05 


5.64 




Hr/M = 


1.25 


5.24 


4.55 


3.93 


20.8 


13.4 


9.57 


8.44 



a ^ ""cent = n cen t/no is the central density ratio of a neutron star containing 
hyperons and nucleons in the absence of kaons. It can be compared with the 
threshold density u c for kaon condensation. 
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